Destruction of small-scale dark matter clumps 
in the hierarchical structures and galaxies 
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A mass function of small-scale dark matter clumps is calculated in the standard cosmological 
scenario with an inflationary-produced primordial fluctuation spectrum and with a hierarchical 
clustering. We take into account the tidal destruction of clumps at early stages of structure formation 
starting from a time of clump detachment from the Universe expansion. Only a small fraction of 
these clumps, ~ 0.1 — 0.5 %, in each logarithmic mass interval AlogAf ~ 1 survives the stage 
of hierarchical clustering. The surviving clumps can be disrupted further in the galaxies by tidal 
interactions with stars. We performed the detailed calculations of the tidal destruction of clumps 
by stars in the Galactic bulge and halo and by the Galactic disk itself. It is shown that the Galactic 
disc provides the dominant contribution to the tidal destruction of small-scale clumps outside the 
bulge. The results obtained are crucial for calculations of the dark matter annihilation signal in the 
Galaxy. 

PACS numbers: 12.60.Jv, 95.35.+d, 98.35. Gi 



I. INTRODUCTION 

One of the promising indirect manifestation of the 
Dark Matter (DM) particles is their possible annihila- 
tion in the Galactic halo 01 . A local annihilation rate is 
proportional to the square of number density of DM par- 
ticles. Therefore, annihilation proceeds more efficiently 
in the dense DM substructures of the Galactic halo. Both 
analytical calculations 2, 3, .4, 5] and numerical simula- 
tions IE 13 with the inflationary-produced adiabatic 
density fluctuations predict the existence of DM clumps 
in the Galactic halo. The enhancement of the annihila- 
tion signal due to the presence of substructures in the 
Galactic halo depends on the fraction of the most dense 
small-scale clumps 13, lo^ . The most essential characteris- 
tics of clumps for calculations of DM annihilation in the 
Galactic halo are the minimum mass and distribution 
function of clumps. At the same time the tidal destruc- 
tion of clumps [j| strongly influences the number density 
of clumps in the Galaxy. 

The small-scale clumps 13, Isl, 1^ are formed only if the 
corresponding density fluctuations are large enough. The 
inflation models predict the power-law primordial fluctu- 
ation spectrum with an power index Up « 1.0. The small- 
scale clumps are formed earlier than the larger ones and 
captured by the larger clumps in the process of a hierar- 
chical clustering in the expanding Universe. Eventually 
all clumps consist in part of the smaller ones and of the 
free DM particles. An effective index of the density per- 
turbation power spectrum ti — > — 3 at small-scales (when 
mass inside the perturbation M —>■ 0). This means that 



a gravitational clustering of small-scale structures pro- 
ceeds very fast. As a result the formation of new clumps 
and their capturing by the larger ones are nearly simul- 
taneous processes. 

A convenient formalism, which describes statistically 
this hierarchical clustering, is the Press-Schechter theory 
[l^ and its extensions, in particular the 'excursion set' 
formalism developed by Bond et al. jl3| (for a clear intro- 
duction see ^3)- However, this theory does not include 
an important process of the tidal destruction of small 
clumps inside the bigger ones. This process has been 
taken into account in our previous work Q , where it was 
demonstrated that only a small fraction of the small-scale 
clumps survives the tidal destruction in the hierarchical 
clustering. Nevertheless, even this small fraction of sur- 
vived small-scale clumps is enough to dominate the DM 
annihilation rate for the most reasonable spectra of pri- 
mordial fluctuations. 

A mass distribution of small-scale clumps survived in 
the hierarchical structuring was derived in Q : 
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where M is a clump mass, n is a power-index of density 
perturbations at a mass-scale M. The distribution func- 
tion ^int is a mass fraction of DM in the form of clumps 
in the logarithmic mass interval dlogM. 

The minimal mass of DM clumps Mmin is determined 
by the leakage of DM particles from the growing den- 
sity fluctuations (the diffuse leakage and free stream- 
ing ) a nd depends on the properties of DM particles 
[a tH El lli El- The existing estimates of 

Afmin for neutralino DM are substantially different, from 
MpitL^ lO-i^A^i in [3 to Mmin ~ (10-"^ - \^-^)Mq in 
[iSl Ha . [l7| . In 01 we performed detailed calculations of 
a DM particle diffusion and free streaming in the kinetic 
equation approach. For the case of neutralino considered 
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as a pure bino we obtained for the minimal mass of DM 
clumps 




where — fr? — m^, with m-^ being a neutralino mass, 
TO is a sfermion mass, and A* has only a logarithmic 
dependence on M and to^. In the considered range of 
parameters A* ~ 83. Our value of Mmin agrees reason- 
ably well with |0, 0, and strongly disagrees with 
hj for the reasons explained in 

Due to uncertainties in the SUSY parameters, a nu- 
merical value of Afniin is not exactly predicted. With our 
choice of SUSY parameters 3, Afmin ~ 1O"^M0, is of 
the Moon-scale mass. With other choice of SUSY param- 
eters 0: Afniin ^ 1O~^M0, is of the Earth-scale mass. 
In the numerical simulations [23, Afmin ~ lO~^Ar0, was 
actually assumed by putting the corresponding cut-off in 
the initial density perturbations. 

The very interesting numerical simulations of the for- 
mation of small-scale DM clumps with a mass larger than 
the Earth mass have been performed recently in |20j |. 
There is direct correspondence of this simulations with 
the earlier theoretical calculations: 

(i) The density profile of large-scale clumps is influ- 
enced by the hierarchical clustering of the smaller ones. 
The new important result of numerical simulations pol | 
is a resolution of the density profile of the isolated min- 
imal mass clumps, Mmin ~ IO^^A/q. The clumps of 
minimal mass are formed directly from the isolated fluc- 
tuations and their density profile is not influenced by the 
hierarchical clustering. The internal density profile of 
small-scale clumps in these simulations is proved out the 
same as in the theoretical calculations |23| performed for 
the isolated density fluctuations. The agreement between 
the theory and numerical simulations for the predicted 
internal density profile of clumps, oc r~^, is fairly good 
within the involved uncertainties: /3 — 1.7— 1.8 in |2ll |. 
and /? = 1.5 — 2.0 in |0 respectively. 

(ii) The numerical simulations |23| agree rather well 
with the shape of theoretically derived mass function of 
small-scale clumps but with the different normaliza- 
tions. 

A tidal destruction in the Galaxy of the Earth-size 
clumps from simulations '2(/| has been recently consid- 
ered in [23. |23- 24. 1 . The results are rather controver- 
sial. Authors of |22ll23| conclude that all the Earth-mass 
clumps are destructed in tidal interactions with stars in 
the Galaxy, while in this result was not confirmed 
under a different assumption on the star number density. 

In this paper we present the alternative and indepen- 
dent calculations for all processes of the tidal destruction 
of small-scale clumps: (i) in the hierarchical clustering, 
(ii) by stars from the stellar bulge, (iii) by stars from the 



halo and (iv) by the Galactic disc. The last of these pro- 
cesses is turned out the most effective. We also describe a 
new method for calculations of clump destruction in the 
hierarchical clustering, which is a more general (valid for 
the arbitrary spectra of primeval fluctuations) and for- 
mally more transparent than the earlier one in 

Our calculations of the tidal destruction of clumps by 
stars in the Galaxy are quite different from Q2sil^2M> by 
both, methods and results. While in the references above, 
only the tidal destructions in collisions of clumps with the 
individual stars were studied, we found that dominant 
effect is provided by the destruction of clumps in the 
collective gravitational field of the Galactic disc. As a 
result we predict that only 17% of the Earth-mass clumps 
survived the tidal destruction at the position of the Sun. 
This result is crucial for the rate of DM annihilation in 
the Galaxy. 

The paper is organized as it follows: In the Sec. ITU we 
describe a new method for calculation of the small-scale 
clump destruction in the hierarchial clustering. We cal- 
culate a mass function of survived clumps and compare 
it with a corresponding one from numerical simulations. 
In the Sec. Illll the tidal destruction of clumps by the 
Galactic disk is considered. In the Sec. II VI the life-time 
of clumps in the central stellar bulge and stellar halo 
spheroid is calculated. In the Sec. we discuss the ob- 
tained results. 

We perform our calculations for the standard cosmo- 
logical model with a matter density parameter — 0.3, 
a cosmological constant term i^A = 1 — f^m — 0.7 and the 
Hubble constant h = 0.7. 



II. DESTRUCTION OF CLUMPS IN 
HIERARCHIAL CLUSTERING 

The process of hierarchical clustering and tidal de- 
struction of DM clumps can be outlined in the follow- 
ing way. The DM clumps of minimal mass are formed 
first in the expanding Universe. The clumps of larger 
mass, which host the smaller ones are formed later, and 
so on. Some part of clumps are destroyed tidally in the 
gravitational field of their host clumps. 

In this Section we study the destruction of DM clumps 
in the process of hierarchical structuring long before the 
final galaxy formation. At small-mass scales the hierar- 
chial clustering is a fast and rather complicated nonlinear 
process. We use a simplified model which nevertheless 
takes into account the most important features of hier- 
archial clustering. 

To describe the formation of clumps we will use the 
model of spherical collapse IT5| in flat cosmology without 
the A-term. This assumption is well justified at early 
times of clumps formation when the A-term is negligible 
in comparison with the matter density. In this model a 
formation time of clump with an internal density p is t = 
{npcj pY^Hcc^, where k = 187r^ and peq = Po(l + -Zoq)^ 
is a cosmological density at the time of matter-radiation 
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equality ^cq, 1 + ^oq = 2.35 x lO^rim^^ and po = 1-9 x 
10~^^nmh^ g cm^'^. The index 'eq' here and throughout 
below refers to quantities at the time of matter-radiation 
equality icq- 

The DM clumps of mass M can be formed from den- 
sity fluctuations of different peak-height v — Scq/(Jcq{M), 
where acq{M) is the fluctuation dispersion on a mass- 
scale M at the time tcq. A mean internal density of 
clump p is fixed at the time of clump formation and 
according to ^3 equals p = Kpcq[i'crcq{M)/Sc]'^ , where 
Sc = 3(127r)2/3/20 ~ 1.686. 

The tidal destruction of clumps is most effective at 
the early epochs of the Galactic halo formation, when 
the host density profiles are not finally established. The 
tidal interaction of clumps is a complicated process and 
depends on many factors: a clump formation history, 
host density profile, an existence of different substruc- 
tures inside the host, orbital parameters of individual 
clumps in the hosts, etc. Only in numerical simulations 
all these factors can be taken into account properly. In 
this paper we use a simplified approach by calculating an 
energy gain per each tidal interaction and a number of 
tidal interactions per dynamical time in the hosts. 

An internal energy of self-gravitating object increases 
in tidal interactions. This energy increase was calculated 
e. g. in for the case of a star globular cluster in a 
spherical galaxy. By using the model of tidal heating 
from _26] , we determine now a survival time (or a time of 
tidal destruction) T of some chosen small-scale clump due 
to the tidal heating inside of a larger mass host clump. 
The motion of a clump would be rather complicated in 
the case of a fast hierarchical clustering of hosts. Dur- 
ing a dynamical time in the host tdyn — 0.5(Gpft,)^^/^, 
where ph is a mean internal density of the host, the cho- 
sen small-scale clump can belong to several successively 
destructed hosts. We will consider a typical clump orbit 
inside the host and assume for simplicity in this Section 
the isothermal internal density profile of the clump. 

A clump trajectory in the host experiences successive 
turns accompanied by the "tidal shocks" . For 

the considered small-scale clump with a mass M and ra- 
dius R, the corresponding internal energy increase after 
a single tidal shock is 



47r 

AEc^—-fiGphMR^, 



(3) 



where a numerical factor 71 ~ 1. Let us denote the 
number of tidal shocks per dynamical time tdyn by 72. 
A corresponding rate of clump internal energy growth is 
E ~ 72Ai?/idyn- A clump is destroyed in the host if 
its internal energy increase due to tidal shocks exceeds a 
total energy \E\ ~ GM'^/2R. As a result, for a typical 
time T = T{p, ph) of the tidal destruction of a small-scale 
clump with density p inside a more massive host with a 
density ph we obtain: 



T-\p,pu)^E/\E\ 



47i72Gi/V'/'p' 



It turns out that a resulting mass function of small-scale 
clumps (see in this Section below) depends rather weakly 
on the value of 7172. 

During the lifetime of an individual small-scale clump, 
it can sequentially inhabit in many host clumps of larger 
mass. After the tidal disruption of the first lightest host, 
a small-scale clump becomes a constituent part of a heav- 
ier one, etc. The process of hierarchical transition of a 
small-scale clump from one host to another occurs al- 
most continuously in time up to the final host formation, 
where the tidal interaction becomes inefficient. 

A corresponding mass fraction of small-scale clumps 
with mass M escaping the tidal destruction in hierarchi- 
cal clustering (or probability of clump survival) is given 
by the exponent function e"'' with 



J: 



E 



T{p,Ph)' 



(5) 



Here Ai/j is a difference of formation times tt of two suc- 
cessive hosts, and summation is over all clumps of inter- 
mediate mass-scales, which successively host the chosen 
small-scale clump of a mass M . By changing the sum- 
mation to integration in (jSJ we obtain 



dth 



T{p, Ph) 



1 



Pi ~ Pf 



Pi 



where 



7 = 27i72«:i/2Gi/Vcq^'ioq ^ 14(7i72/3), 



(6) 



(7) 



and ti , tf, pi and p / are respectively the formation times 
and internal densities of the first and final hosts. It can 
be seen from ^ that the first host provides a major con- 
tribution to the tidal destruction of small-scale clumps, 
especially if the first host density pi is close to p, and 
consequently e^'^ ^ 1. 

A mass function of small-scale clumps (i. e. a differen- 
tial mass fraction of DM in the form of clumps survived 
in hierarchical clustering) can be expressed as 



i^diy = dMdv{2/Ti)e-''''/'^ 



dvi e 



X I dt 



d^F{M, i) 



DM dt 



^Jlp(u)Xi)] 



(8) 



(4) 



In this expression to is the present age of the Uni- 
verse and F{M, t) is the mass fraction of unconfined 
clumps (i. e. clumps not belonging to the more massive 
hosts) with a mass smaller than M at time t. Accord- 
ing to ^2 1 the mass fraction of unconfined clumps is 
F{M,t) = eTi{SJ[V2aeq{M)D{t)]), where erf(a;) is the 
error-function and D{t) is the growth factor normalized 
by D(tcq) = 1. An upper integration limit to in is 
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not crucial and may be extrapolated to infinity because 
a main contribution to the tidal destruction of clumps 
is provided by the early formed hosts at first steps of 
hierarchical clustering. 

Two processes respond for a time evolution of the frac- 
tion dF{M, t) I dM dM of unconfincd clumps in the mass 
interval dM: (i) the formation of new clumps and (ii) 
the capture of smaller clumps into the larger ones. Both 
these processes are equally effective at the time when 
d'^F /{dMdt) = 0. To take into account the confined 
clumps (i. e. clumps in the hosts) we need only the 2nd 
process (ii) for the fraction dF(M,t)/dM . Nevertheless, 
in ^ it is used the fraction dF{M,t)/dM influenced 
by the both processes. This is not accurate at a typical 
formation time of clump with a mass M, when clump 
density is comparable with the density of hosts. Fortu- 
nately, for this time the exponent in ((SJ is very small, 
e~'^ <C 1, as it can be seen from (0 and iQ. Respec- 
tively, an uncertain contribution from the process (i) to 
the integral ||HJ) is also very small. Meanwhile, only the 
process (ii) dominates in the integration region where the 
exponent e~'^ is not small. For this reason ||SJ| provides 
a suitable approximation for the mass fraction of clumps 
survived in the hierarchical clustering. 

Finally we transform the distribution function © to 
the following form: 



^ M 



2Tr ' dM 



(9) 



Here the numerically calculated function y{v) depends 
rather weakly on the parameter 7 from Q and is shown 
in the Fig.^ By deriving lO, we take into account that 
a{M) is a slowly varying function of clump mass M . For 
the same reason, providing an integration in ijSJ we use 
the dependance of ti(pi) only on the variable v by ne- 
glecting the dependance on M . Physically the rising of 
y{i') with V corresponds to a more effective survival of 
high-density clumps (i. e. with large values of v) with 
respect to the low-density ones (with small values of v). 
Integrating ^ over we obtain 



^ dM , , dM 

6„t^-0.017(n + 3)— . 



(10) 



This mass function is in a reasonable agreement with 
the similar one from in the our earlier work An 
effective power-law index n in IjlUfl is given by 



1 



, c)l0gCrcq(M) 

d log M 



(11) 



and depends very weakly on M . At the small mass-scales 
one has approximately n ~ rip — 4. Equation implies 
that for the suitable values of n only a small fraction of 
clumps, about 0.1 — 0.5 %, survives the stage of hierarchi- 
cal tidal destruction in the each logarithmic mass interval 
AlogM ~ 1. A simple M~^ shape of the mass function 
(|10|l is in a very good agreement with the corresponding 
one obtained recently in the numerical simulations pol | . 



> 



FIG. 1: Numerically calculated function y{v) from 



It must be stressed that a physical meaning of the sur- 
vived clump distribution function ^ dv dM /M is different 
from the similar one for the unconfined clumps, given 
by the Press-Schechter mass function dF/dM . For com- 
parison, the Press-Schechter mass function of unconfined 
clumps 12] is 



eps(i) 



dM 



da, 



cq 



2TTa^^D{t) dM 



exp 



2<^^(i) 



dM, 



(12) 

where a^q — foq(Af) The mass function of clumps sur- 
vived in a hierarchical clustering (|10|l is several times less 
than the Press-Schechter mass function H12|) at a mean 
time of clump formation with acq{M)D{t) ~ Sc- 

In further calculations we will use an interpolation fit- 
ting of the fluctuation dispersion CTeq(M) from (27|| (see 
also 



(M) 



2 X 10-* 



k 

km 



(np-l)/2 



"-cq 



1 3/2 



(13) 

where the wave vector k oc M *''^, respectively k^q and 
km correspond to a mass inside the cosmological horizon 
at the moments icq and to, Up is a primordial pertur- 
bation index, and fsi^K) = 1-04 - 0.82f]A + 2f7^. It 
must be noted that interpolation (|13() is valid only for 
the small-scale clumps, with M < 10'^ Mq. The analy- 
sis of the WMAP data of the CMB anisotropy re- 
veals a power-law spectrum of initial perturbations with 
Hp — 0.99 ± 0.04 in a good agreement with the canon- 
ical inflation value Hp — 1.0. However, when the data 
from 2dF galaxy power-spectrum and hy-a are included 
in the analysis, the best-fit favors in a softer spectrum 
with Up = 0.96 ± 0.02. Nevertheless, the recent observa- 
tions does not exclude even the values Up — 1.1. 

Note that a differential number density of small-scale 
clumps in the Galactic halo n{M) dM cx dM/M^ from 
(|10|l coincides not only with a similar one from the re- 
cent numerical simulations of small-scale clumps |23| but 
is also very close to that obtained in the numerical sim- 
ulations for large-scale clumps with mass M > IO^Mq. 
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See Fig. [5]for a comparison. 

Strictly speaking, our calculations are not valid for 
large-scale clumps because of their continuing tidal de- 
struction in the halo up to the present epoch to and the 
accretion of the additional large-scale clumps into the 
halo from the intergalactic space. Nevertheless, our ap- 
proach remains valid even for the large-scale clumps in 
the narrow mass range, where the power-law perturba- 
tion spectrum can be used as a rather good approxima- 
tion. 

In the Fig. [21 a differential number density of small- 
scale clumps from (I10|) is shown by the solid line. As 
it was noted above, the region of validity for this curve 
is M < 10'^ Mq. For larger masses an extrapolation is 
shown (right part of the solid line). The corresponding 
mass functions from numerical simulations can be param- 
eterized in the form £_{M)dM/M = AM^-^dM. The 
constant A can be determined by fixing a power-law index 
A and a fraction e of the halo mass in the form of clumps 



with mass from A/n 



M„ 



£ = 



10'' M© to M„ 



AM^-^ dM. 



10i°Mp 



©• 



(14) 



M„ 



With this parametrization a mass function of large-scale 
DM clumps may be expressed as 



C(M) 



dM 



dM 
IF 



(2 - \)M 



2-A 



2-A 



mm 



J^-' mm 



A ^2; 



A = 2. 



(15) 



In the Fig. [3 a differential number density of large-scale 
clumps (A/h/M) i{M) dM/M from ^ is shown for dif- 
ferent values of A and e taken from various numerical 
simulations: e ~ 0.2, A = 2 from 0; e ~ 0.15, A = 1.9 
from 0, and e ~ 0.05, A = 1.78 from [s^. Obser- 
vations of the Galactic halo lensing |^ give a smaller 
clump fraction value, e ~ 0.02. One can see in the Fig. [21 
a reasonable agreement between the extrapolation of our 
calculations and the corresponding numerical simulations 
of the large-scale clumps. 



III. DESTRUCTION OF CLUMPS BY DISK 

Crossing the Galactic disc, a clump can be tidally de- 
structed by the collective gravitational field of stars in the 
disc. This phenomenon is similar to the destruction of a 
globular cluster by the "tidal shocking" in the Galactic 
disc m. 

The rate of energy gain per unit mass due to the tidal 
shocking was calculated in : 



dE 
H 



4g,l(A^) 



(16) 



1020 



,15 



^ 1 n 



£=0.2, A=2 
e=0.15, A=1.9 
e=0.05, A=1.78 

Equation (10) 



10" 
M/Mo 



,10 



FIG. 2: A differential number of small clumps in the Galaxy 
from Hl()|l for rip = l.Q is shown by the solid line. These calcu- 
lations are valid only for small-scale clumps with M < 10'' Mq 
and extrapolated to the larger masses (the right part of the 
solid line). Other curves are the corresponding number den- 
sities of large-scale clumps with M > from numerical 
simulations (see details in the text) for different values of pa- 
rameters £ and A from 11511 . The left parts of these curves are 
the extrapolations to small masses. 



Here gm is the maximum gravitational acceleration ac- 
quired by the constituent DM particle of the clump mov- 
ing in the gravitational field of the disk, Az = is a vertical 
(perpendicular to the disk plane) distance of a DM par- 
ticle from the clump center, Tc is an orbital period of 
clump in the halo, Vz^c is a vertical velocity of disk cross- 
ing. In 1(1^ the two crossings of disk by a globular cluster 
during the orbital period Tc is assumed, while only one 
disk crossing is typical in the case of elongated orbits of 
DM clumps in the halo (see below). 

A surface mass of the Galactic disk can be approxi- 
mated by a simple exponential law [s^ 



Md 



-r/ro 



(17) 



with Md = 8 X lO^"Af0 and tq = 4.5 kpc. The maximum 
gravitational acceleration during the disk crossing is 



(18) 



Following to |21| , we use a power-law parametrization of 
the internal density profile of DM clumps: 



3-/3 /rx-/5 
Pint(r) = ^p(-j 



where p and R is respectively a mean internal density and 
radius of clump, P — 1.7 — 1.8 and we put Pint('') = 
at r > i?. The corresponding power-law profile of small- 
scale clumps with (3 ~ 1.5 — 2 has been recently found 
in numerical simulations [23|. For this profile a total 
(kinetic and potential) energy of a clump is given by 



\E\ 



3-P GM^ 
2(5 - 2(3) R 



(20) 
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Integrating (|16|l over the clump volume with the profile 
(|19|l . one can obtain a total rate of energy gain by clump 
dE/dt and then a time of clump destruction by the disk: 



M 

E 



2(5-/3) GT.pv 



2 

z,c 



3(5 - 2/3) g2 



(21) 



Note that the adiabatic correction for the disk shocking 
(see e. g. [s^) is very small in the case of DM clumps 
and may be neglected. 

To estimate the tidal shocking effect produced by dif- 
ferent parts of the Galactic disk at radial distance r, let 
us consider at first a toy halo model by assuming the cir- 
cular orbits of DM clumps. Then a disk crossing velocity 
Vz,c equals to a circular velocity: 



Vrot{r) 



1/2 



(22) 



We suppose here for simplicity a pure isothermal den- 
sity profile of the halo 



Ph(0 



4tt Gr^ ' 



(26) 



where wh = (GAfn/^n)^^^ is the halo rotational velocity 
and i?H is the Galactic halo radius. At r > i?H we put 
P'ii{r) = 0. In the absence of an analytical model of the 
finite isothermal sphere, we construct a simplified model 
which approximates the isothermal sphere in the inner 
region, at r ^ i?H- A gravitational potential energy 
U{r) corresponding to the density profile H26|l is 



U{r)=mvl[\og[T/R^)-l], 



(27) 



where m is a mass of DM particle. The radial motion of 
a particle with mass m and angular momentum L in the 
spherical potential obeys the equation 



where M-H_{r) is a halo mass inside the sphere of radius r. 
Using for a clump orbital period Tc — one 
finally finds 



(5-/3) r^pM^/^iry/^ 



37r(5- 2/3) 



G1/2M2 



,2r/ro 



(23) 



Comparison of a clump destruction time from l|23|) 
with the Universe age to shows that all clumps with the 
internal density p < 2 x 10~^^ g cm~^ are effectively 
destructed within the radius r < 15 kpc from the Galactic 
center. In particular, the major part of the Moon-mass 
clumps with M = 2 x 10~^Mq, Up = 1 and v = 2 do not 
survive inside the central 15 kpc. 

However, in the real Galactic halo the DM clumps have 
an elongated orbits in general. These orbits cross the 
stellar Galactic disk only once during the orbital period. 
Therefore, we must introduce a factor 2 in (|21|l for these 
orbits. At the same time it is a much more important that 
clumps with elongated orbits have the longer orbital pe- 
riods Tc than in the previous toy model. This elongation 
of orbits significantly increases the probability of clump 
survival. 

As an example let us consider the Galactic halo model 
with an isotropic velocity distribution. This model is ap- 
propriate for the halo formed by the hierarchial clustering 
of clumps. In this model according to (36] the energy dis- 
tribution function f{E) of DM particles is related with 
the density profile of the halo p}i{r) as 



PH(r) = 2^/2^ J ^E - U{r) f{E) dE, (24) 

U(r) 

1 d 



f{E) 



dr dpii{r) 
2y^7r^ dE J - U{r) dr ' 

T{E) 



(25) 



where U{r) is a gravitational potential energy and func- 
tion r = r{E) is defined by the equation U[r{E)] = E. 



By introducing the dimensionless variables 
_ r E I? 



mvt 



y 



(28) 



(29) 



the equation for the turning points, r'^ 
tial H27|l can be written as 



4-2(.- 



logs -f 1). 



0, in the poten- 



(30) 



The derivatives of the left and right sides of this equation 
are equal at s = y^^^- Respectively, the roots s„iinix,y) 
and Sniax(a;, y) of (|30|) satisfy the condition Sniin{x,y) < 
y^^^ < Smax(a;,y). A condition for the existence of the 
solution of (|30|l is a; > (logy — l)/2. The equality in this 
condition corresponds to the circular orbit with Smin = 
Smax- From H28|l one can determine the orbital period: 



Sma; 

Tc{x,y) = 2^ f 

VH J 



ds' 



^2{x- logs' + \)-y/s'- 



(31) 



In the following we will solve H3U|) and find the orbital pe- 
riod Tc{x,y) from (|31|l numerically. Denoting p = cos 6, 
where 9 is an angle between the radius- vector r and the 
particle velocity v, we have 



y = 2(1 -p2)s2(a; - logs + 1). 



(32) 



We find from H25I) the distribution function of particles 
with an energy a; < — 1 by using the density profile H26f) 
with a cutoff at r = Rn: 



1 



,1/2 



25/27r3e Gm^/^Rl 



Fix), 



where 

F{x) = \/2^e~2xg^f 



v/-2(.x+l) 



(33) 



(34) 
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Note that the isotropic distribution function (|33|) repro- 
duces the density profile (|26|l only in the inner halo re- 
gion, at r ^ _Rh- The assumed isotropy of particle dis- 
tribution (i. e. an independence of distribution function 
on the particle angular momentum L) is violated near 
the boundary of the halo, at r ~ i?H- Nevertheless, the 
distribution function (|33|l is adequate for our purpose be- 
cause a tidal destruction of clumps by the Galactic disk 
takes place only in the inner halo region, at r ^ _Rh- At 
the same time, in the considered model the clumps on the 
outer orbits, at r ~ i?H, provide only small contribution 
to the halo density at r <C i?H- Neglecting these outer 
clumps, we can define the probability of clump survival 
in tidal destruction by the Galactic disk (or the fraction 
of survived clumps) as a function of radius r = sRn in 
the following form: 

1 -1 

J dp J dx - log s + 1 F{x) e"*"/*'' 

P,{r) = ° ■ (35) 

/ dx \Jx — log s -I- 1 F{x) 

log s — 1 

Here td is from l|21|) but with an additional factor 2 
(one disk crossing per orbital period) and with the re- 
placements Tc ^ Tc{x,y) from (EU, .g„i ffm(fmin), 
Vz,c v{r^in), where r,nin = -Smin-^H, gmin I S the mini- 
mal root of (1201) and v{r) = y/2[E - U{r)]/m. See in the 
Figs. 01 |S1 the resulting probabilities of clump survival in 
the Galaxy (or fractions of clumps survived the tidal de- 
struction) . 

Note that the probability of clump survival Pdir) in 
(|35|l is an approximate expression averaged over an an- 
gle between the plane of the Galactic disk and the clump 
orbit plane. In the real Galactic halo there must be some 
anisotropy in clump distribution with respect to the disk 
plane. For example, the clumps with orbits in the Galac- 
tic disk plane are destructed more efficiently than ones 
outside the Galactic plane. 



IV. DESTRUCTION OF CLUMPS BY STARS 

The internal energy increase of a clump during a single 
star flyby is 



AS= d^r pint{r)ivz - Vzf 



(36) 



where Vz is a velocity increase of constituent DM particle 
inside a clump in the direction of axis z and Vz is a similar 
one for a clump center-of-mass. The axis z is directed 
along the line connecting a clump center-of-mass with 
a star at the moment of a star closest approach. In the 
impulse approximation, by neglecting the internal motion 
of DM particles in a clump during the star encounter and 
assuming the straight line orbit of a star (see e. g. p^ ) 
we have 

- A/ = ?^rcosV, (37) 




M/Mo 



FIG. 3: The fraction h/R according to 1421 1 as a function 
of clump mass M at the distance from the Galactic center 
r = 2 kpc, 1/ — 2 and rip = 1.0 and Hp — 1.1 respectively. 



where / is the distance of a star closest approach to a DM 
clump and "0 is a polar angle in the spherical coordinates. 

Let us Wrei is relative velocity of a star with respect 
to a DM clump. In the approximation of a rectilinear 
motion, an angle between the line connecting a clump 
center-of-mass and Wvpi evolves as 



dt 



- —j- COS ( 



(38) 



Changing a variable i to in the Newton equation of 
motion, one gets 



dv. 



Gr 



COS( 



(39) 



where to* is a typical star mass. After integration of this 
equation we obtain 



2G'to* 



(40) 



Now by integrating H36|l over a clump volume with a den- 
sity proffie Pint{r) from H19|l . we find in the case oil > R: 



AE 



2(3-/3) G'^MR'^ 
3(5 - /3) 



w2 J4 

rcl 



(41) 



The opposite case I < R was considered e. g. in It is 
easily verified that the maximum internal energy increase 
occurs for a star flyby with I ~ R. 

At this step we must distinguish two cases: (i) clump 
destruction during a single star flyby and (ii) clump de- 
struction after numerous star collisions. In the first case a 
threshold for clump destruction is achieved at AE = \E\, 
where a total energy of clump E is given by (|20|l . From 
the equality AE = \E\ one finds the maximal impact 
parameter for a single flyby destruction 



dl 



dl 



4(5-2/3) Gml 
3(5-/3) MRvl 



(42) 
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FIG. 4: The fraction of clumps with mass M = 2 x lO'^M© 
and peak-height v = 2 survived a tidal destruction in the 
Galactic disc Pd, in the Galactic halo Ph and the resulting 
total fraction Ptot = P^iPd as a function of distance from the 
Galactic center. The cutoff at r < 3 kpc is due to destruction 
of clumps inside the bulge. 



where V ~ [GM j R)^!"^ is a velocity dispersion of DM 
particles in the clump. The fraction l^jK as a function 
of clump mass M is shown in the Fig. |21 Note that 
condition > 1 is satisfied for clumps of the smallest 
mass. A total rate of clump destruction by stars in the 
case of > 1 is given by 

.-1 ^ ,2 E{l>h) 

= ^ = <n,v,,i + -^J^- (43) 

where n* is a number density of stars and 



FIG. 5: The same as the Fig.Hbut for M = 10"'' Me 
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FIG. 6: The same as the Fig.Hbut for M = 10"^M(r 



E{1 > U) = 27r / A£:(/) n^Vrci I dl. 



(44) 



After integration in H44() with U from (|42() . we find that 
the second term in (|43|) is equal to the first one. Thus, 
the resulting time of clump destruction in the case of 
h/R > 1 is 



1 



3(5-/3) M 
(5 - 2(3) GR^ 



1/2 



(45) 



We see from (I45f) that the time of clump destruction by 
stars does not depend on v^d in the case of > 1. 

Similarly, the time of clump destruction by stars in the 
case of U < R is 



3(5 - /3) v,,iM 
87r(5 - 2/3) GRmln 



(46) 



A. Destruction of clumps in the bulge 

The bulge is central spheroidal subsystem of the 
Galaxy. Following to |^ we approximate the radial 



number density distribution of stars in the bulge in the 
radial range r = 1 — 3 kpc as 



nbA'^) = ipb/m^)e^p [-{r/n) 



1.61 



(47) 



where pb — BMq/pc'^ and rb — I kpc. By using H47(l 
together with (|45|l or (|46|l it can be shown that inside 
the bulge, at r < 3 kpc, all small-scale clumps with M > 
10~^Mq are tidally destructed during the Hubble time, 
i. e. <C tQ. Therefore, there is an empty cavity in 
clump distribution in the Galactic center with a radius 
r ~ 3 kpc as it can be shown in the Figs. 

What is a fate of the core of a tidally destructed 
clump? Let us consider the scaling of destruction time t^, 
in dependence on a varying clump radius r and mass 
M{r) cx r^^^. According to H45|l or (|46|l . a time of 
clump destruction is scaled respectively as (x r~^/^ 
or (X r^~^. For surviving of the core, it is necessary 
that ^ at r "> 0. This is possible only if the in- 
ternal density profile p{r) (x is rather steep, /3 > 2. 
Meanwhile, both the theoretical models and numerical 
simulations predict /3 < 2, and therefore, the core does 
not survive during the tidal destruction of DM clump. 
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B. Destruction of clumps in the halo 

The radial number density distribution of stars in the 
Galactic halo (outside the Galactic disk) at radii r > 
3 kpc can be approximated as 

rihAr) = ipQ/m,){rQ/rf, (48) 

where p© = 10~^ Mq /pc^ and = 8.5 kpc. The stel- 
lar density profile in the Galactic halo is rather poorly 
known, and so H48(l must be considered only as an upper 
limit (2J|. We will describe the distribution of clumps 
in the halo the same way as in the Section ITTll To take 
into account a varying number density of stars n/i.*(r), 
we made an averaging of the rate of clump destruction by 
stars t^^ along the orbital trajectory during the orbital 
period: 

(t-if^ )]^ T ds't:^ 

U {x,y)) y^T,ix,y)J ^2{x -logs' + l)~y/s'^' 

(49) 

In this expression the dimensionless variables x and y are 
from l|29(l . a clump orbital period Tc{x, y) is from (|31() and 
destruction time of clump is from (|45|) or H4()|) with the 
replacements Us ^ nh^*{r) and =^ \/2[i? — U {r)]/m. 
The averaging procedure H49|l is in fact the integration of 
energy gain rate J Edt along the clump orbit. 

The resulting probability of clump survival in tidal 
destruction by the Galactic halo stars P}i{r) is defined 
by a similar expression as H35|l but with a replacement 
g-to/td ^ e-to{c'(^,y)>^ where {t-^x,y)) is from 

The results of numerical calculations of the tidal de- 
struction of DM clumps by different Galactic components 
are summarized in the Figs.EHEl These calculations were 
performed for DM clumps originated from fluctuations 
with the peak- height v ~ 2. 

Correspondingly 32%, 27% and 18% of clumps survive 
the destruction by the Galactic disk tidal shocking at 
the Sun position, r© — 8.5 kpc, for clump masses M = 
2 X 10-®Mq, M = IO-^Mq and M = IQ-^Mq. The 
Galactic disk destroys clumps even outside its bound- 
ary, at r > 15 kpc, because some of DM clumps with 
the extended orbits intersect the Galactic disk in the in- 
ner part of the halo. The destruction of clumps by the 
Galactic disk becomes inefficient at r > 40 kpc. The re- 
spective fractions of clumps of the same masses surviving 
the tidal destruction by stars in the Galactic halo (out- 
side the Galactic disk) are 66%, 63% and 57%. The final 
fractions of clumps of the same masses survived the tidal 
destruction both by the Galactic disk and stars in the 
Galactic halo F(rQ) = PH(ro)Pd(ro) are 21%, 17% and 
10% respectively. 

In the Fig.[7|the fraction of survived small-scale clumps 
in the Galactic halo is shown in dependance on a mean 
internal density of small-scale clumps. 




FIG. 7: The survived fraction of small-scale clumps Ptot{r) in 
the Galactic halo inside the radial distance r < 100 kpc. A 
mean internal density of clump p is in g cm"''. 



V. CONCLUSIONS 

We calculated the number density distribution of 
small-scale DM clumps in the Galactic halo in depen- 
dance on a clump mass M, radius R (expressed through 
the fluctuation peak-height v) and radial distance r to 
the Galactic center. These calculations were performed 
by taking into account the tidal destruction of clumps in 
the early hierarchical clustering and later in the Galaxy. 

Calculations of the distribution function of small-scale 
clumps are carried out, following to 4], in the framework 
of the standard cosmological model and the hierarchical 
model of structure formation. The primeval power spec- 
trum of density perturbations P{k) (x is taken from 
the inflation models with rip ~ 1 (the Harrison- Zeldovich 
spectrum). In this model the small-scale clumps are 
formed earlier than the bigger ones. The minimal mass 
Mmin of clumps is determined by the free streaming of 
DM particles from a growing fluctuation. The value of 
Mmin is a model dependent quantity. For neutralino as 
DM particle, the minimal mass Mmin is given by lj2Jl, and 
it is the Moon-scale mass. 

In the process of hierarchical clustering the small 
clumps are captured by the bigger ones, and so on. Thus 
the hierarchical structure is formed, when all clumps con- 
sist in part the smaller ones and the free DM particles. 
Some part of DM clumps are tidally disrupted in the 
gravitational field of the bigger host clumps. In this sce- 
nario we calculated the differential distribution of the 
survived clumps given by |(5J) as a function of two inde- 
pendent parameters: e. g. a clump mass M and fluc- 
tuation peak-height v (or a clump mass M and radius 
R). The corresponding integral mass function is given 
by (fTn|l . where the small factor ^int — 0.017(n -I- 3) gives 
the mass fraction of clumps survived the tidal destruction 
in the hierarchical structuring. 

The predicted differential number density of small 
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clumps £,{M)(j)ji/M) dM/M is very close to our previous 
calculations 0], and both are in a good agreement with 
the recent results of numerical simulations [2(tI |. 

Our calculations are valid only for small-scale clumps 
with masses M < 10"^ Mq. The physics of larger mass 
clumps is rather different. For large-scale clumps the 
dynamical friction, tidal stripping and accretion of new 
clumps into the halo proceed in a different way. Never- 
theless, the calculated mass function is in a good agree- 
ment with a mass function of the large clumps (obtained 
in the numerical simulations) in the intermediate mass 
range (see in the Fig.|21). 

The mutual tidal destruction of small-scale DM clumps 
is effective only at the early stage of hierarchical cluster- 
ing. At later stages the DM clumps are additionally de- 
structed by stars and by the collective gravitational field 
of the Galactic disc. In the Galaxy at radial distance 
r < 3 kpc all small-scale clumps are destructed by stars 
in the central bulge. At radial distances in the range 
r = 3 — 40 kpc the DM clumps are destructed by stars 
from the halo and by the tidal shocking in the Galactic 
disk. The latter provides the major contribution to the 
tidal destruction of clumps outside the bulge. Only 21%, 
17% and 10% of clumps survive the tidal destruction near 
the Sun position for clump masses M — 2 x IO^^A/q, 
M = 10~^Mq and M = IO'^Mq respectively Our re- 
sults on the tidal destruction of clumps differ from both 
[23. l23| and (2^, with the intermediate conclusions. At 



radial distances r > 40 kpc the destruction of clumps 
by the Galactic disk becomes inefficient, and the number 
density of clumps is determined only by the early epoch 
of hierarchical clustering. 

The tidal destruction of clumps by the Galactic disk 
and stars affects the annihilating signal mainly in the cen- 
tral region of the Galaxy where destructions are most ef- 
fective. Therefore, a growing fraction of survived clumps 
P(r) smooths the anisotropy of the awaited annihilation 
signal at the Sun position. A local annihilation rate is 
proportional to the clumps number density and, respec- 
tively, to P{r). For example, at the position of the Sun 
the 17% of clumps survive, and so the local annihilation 
rate more then 5 times less in comparison with the P ^ 1 
case. 
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